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Kinetics of water flow through polymer gel 
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Abstract. The water flow through the poly(acrylamide) gel under a constant water pressure is measured 
by newly designed apparatus. The time evolution of the water flow in the gel, is calculated based on 
the collective diffusion model of the polymer network coupled with the friction between the polymer 
network and the water. The friction coefficient are determined from the equilibrium velocity of water 
flow. The Young modulus and the Poisson's ratio of the rod shape gels are measured by the uni-axial 
elongation experiments, which determine the longitudinal modulus independently from the water flow 
experiments. With the values of the longitudinal modulus and of the friction determined by the experiments, 
the calculated results are compared with the time evolution of the flow experiments. We flnd that the time 
evolution of the water flow is well described by a single characteristic relaxation time predicted by the 
collective diffusion model coupled with the water friction. 

PACS. 83.10.Bb Kinetics of deformation and flow - 83.80.Kn Physical gels and microgels 



1 Introduction 

The gel is an important state of matter that is found in 
a wide variety of biological, chemical, and food systems 
pp. The polymer gel consists of a cross-linked polymer 
network and a large amount of solvent (typically water) . 
Since the average mesh size (the distance between neigh- 
boring crosslinks) of the polymer network is in general 
large compared to small molecules, they can pass through 
the gel easily. Two different transport processes are im- 
portant in the gel. 

The first one is the diffusive flow of molecules |2j . The 
transport of molecules by the diffusion in the polymer net- 
work is modeled as the diffusion of the probe objects in the 
fixed mesh of obstacles. The ratio between the probe size 
and the mesh size plays essential role to determine the dif- 
fusion coefficient of the molecules in the gels. The diffusion 
of the probe molecules in the gel has been measured by 
the pulsed field gradient nuclear magnetic resonance. The 
results indicate that the diffusion coefficient of the probe 
molecules of various sizes in the gel is well described by a 
simple scaling relationship [31. 

The second process is the convective flow of the wa- 
ter through the polymer network [HElll] . When the water 
flows in the gel, it experiences the hydrodynamic friction 
from the polymer network, at the same time, the polymer 
network of gel is deformed from the initial conflguration 
by the drag force of the water. The entire flowing pro- 
cess of the water is determined by the balance between 
the viscoelastic response of the water flow and that of the 



polymer network. The collective diffusion of the polymer 
network coupled with water friction, therefore, plays es- 
sential roles in the convective flow process in the gel. 

As well as the transport phenomena in the gel, the 
water flow is also important in the kinetics of the vol- 
ume phase transition of the gel including the pattern for- 
mations in shrinking gels 7,8,9 i[Tni[IIl[Tl[Tl[Tl[Til[T51[T7]. 
Kinetics of solvent flow in the gel is the main subject of 
the application of the gel to control the solvent flow [HI 
[T^ . Although the significance, the solvent flow process in 
the polymer gel has yet to be studied in detail. 

It is the purpose of this paper to describe the kinetics 
of the water flow through the polymer gel, based on the 
collective diffusion model coupled with water friction. The 
characteristic relaxation time of the water flow process in 
the gel is expressed as r = L^/tt^Dc with the parameters: 
the typical size of the gel, L, and the collective diffusion 
coefficient of the polymer network, Dc- The collective dif- 
fusion coefficient is expressed by using the longitudinal 
modulus of the gel, k, and the friction coefficient between 
the polymer network and the water, /, as [511^ 



Dr 



/■ 



(1) 



The water flow in the polymer network is, thus, directly 
governed by the viscoelastic property of the gel. The col- 
lective diffusion coefficient of polymer gels has been ex- 
tensively measured by quasi-elastic light scattering [^ 
I23j . The friction coefficient between the polymer network 
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water gel 




Fig. 1. Experimental setting of water flow measurements 



of gel and the water has been measured by flowing the 
water through the gel [lll5l[6]. 

In this paper, we present experimental results on the 
water flow through the gel, and show that the time evolu- 
tion of the water flow through the gel is described satisfac- 
torily by the collective diffusion model of polymer network 
with water friction. The information on the kinetics of wa- 
ter flow in the gel will help the better use of them in the 
applications. 




Fig. 2. Apparatus for water flow measurements: (a) the up- 
per photograph shows the external appearance, (b) the lower 
photographs shows the inside of the end plates. 



2 Experiment 

2.1 Water flow measurements 

The experimental setting of the water flow measurements 
is illustrated schematically in Figure 1. A thin circular 
slab gel of thickness, L, is set in a tube and the rim of the 
gel is glued to the tube. The surface of the gel at the lower 
reaches side is mechanically fixed to the end plate through 
an 0-ring. When the small pressure, Pq, is applied to the 
water from the upper stream side of the gel, the water 
flows through the gel deforming the polymer network of 
gel by the drag force. The macroscopic value of the friction 
coefficient of the gel, /, is determined from the velocity of 
the water flows out of the gel in the equilibrium state, Uoo, 
as 

Even though the concept is simple, the mechanical 
measurement of the friction coefficient between the poly- 
mer network of gel and the water is not easy because the 
friction of the gel is huge. The difficulty of evaluating the 
friction coefficient measured from the water flow through 
the gel is described in the reference |4]. In our experiments, 
the difficulty is overcome by using a well calibrated homo- 
geneous glass capillary to amplify the velocity of water as 
shown in Figure 1 . The position of the meniscus of water 
in the capillary is measured as a function of time after the 
pressure is applied to the water. 



The photographs of the apparatus used in this study 
are given in Figure 2. The essential structure of the appa- 
ratus is the same with the previous one but some points 
are improved for the precise measurements of the friction 
of the gel 4 . The right hand side of Figure 2 (a) is the 
upper stream side and is connected to a water column to 
apply the hydrostatic pressure. The glass capillary is set 
at the top of stainless steel pipe which can be seen at the 
left hand side of Figure 2 (a) . Leak of water around the gel 
is prevented by setting 0-rings inside the cell as shown in 
Figure 2 (b) . All parts of the cell is made of stainless steel 
to avoid the mechanical deformation due to the applied 
pressure. The mechanical deformation of the apparatus is 
considerably reduced compared to the previous one. 



2.2 Mechanical response measurements 

The mechanical response of the polymer gel relates the 
friction coefficient to the collective diffusion coefficient of 
the polymer network as Dc — k/ f where k = K -\ — ^ 
[20^. Here, k is the longitudinal modulus of the gel and K 
and fi are the bulk modulus and the shear modulus of the 
gel, respectively. We measure two elastic moduli of the gel 
simultaneously, and the longitudinal modulus of the gel is 
determined from them. 

The rod shape gel is elongated in water at a strain of 
about 10 % and the stress is measured. The size of the 
gel is observed by using a microscope during the elasticity 
measurements to obtain the Young's modulus of the gel. 
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E, and the Poisson's ratio, a. The longitudinal modulus 
of the gel, K, is determined from the Young's modulus and 
the Poisson's ratio by k^K+^,K^ E/3{1 - 2a) and 
IJ = E/2{l + <j). 

In this short time experiments, the water solvent nei- 
ther goes out nor comes into the gels. Thus we determine 
the mechanical response of the gel independently from the 
friction experiments by the water flow. We use the results 
to analyze the kinetics of water flow through the polymer 
network. 



2.3 Sample 

The samples used in this paper are poly(acrylamide) gel 
that is obtained by the radical co-polymerization of the 
main-chain component, acrylamide, and the cross-linker, 
N,N'-methylene-bis-acrylamide. The concentration of the 
gel used in this measurement is 1 M. Ammoniumpersul- 
fate and N,N,N',N'-tetramethyl-ethylenediamine are used 
as the initiator and the accelerator. All chemicals used 
here are of electrophoresis grade, purchased from BioRad 
and used without further puriflcation. The gel for the fric- 
tion measurements is prepared in the circular gel mold, 
where the gel-bond films (FMC) with circular opening are 
glued on the both sides. For the friction coefficient mea- 
surements, the gel is prepared in the measurement cell 
and kept under water at a constant room temperature for 
overnight to reach the equilibrium state. 

For the elastic modulus measurements, the rod shaped 
gels are prepared in the capillary. The desired amounts 
of the main-chain component, the cross-linker, and the 
accelerator are dissolved into the distilled and de-ionised 
water, which is prepared by a Mifi-Q system. The pre- 
gel solution is de-gassed for 20 min, and then the desired 
amount of the initiator is added to the solution to initiate 
the reaction. The sample gels are taken out of the reaction 
bath and then washed by distilled and de-ionized water 
extensively and used in each measurement. 



3 Experimental results 

3.1 Friction coefficient of gel 

The experimental results of the measurement of water flow 
through the gel are shown in Figure 3. The position of the 
meniscus in the capillary of lOjil glass micro-pipette is 
plotted as a function of the time elapsed after the appli- 
cation of the pressure to the water. The slope of the flow 
curve, thus obtained, represents the water velocity in the 
capillary, v'. The velocity of the flow is high at first. Then 
it eventually decreases with time and approaches to an 
equilibrium value. The equilibrium state is attained about 
5 to 7 X 10'^ s after the pressure is applied to the water. 
The velocity of the meniscus in the equilibrium state is de- 
termined by the least-squares analysis of the linear portion 
of the flow curve in Figure 3, typically in the time region 
more than 1 x 10"* s, that yields to v'^ = 1.9 x 10^^ m/s. 
The velocity of water flow in the gel, v, is then calculated 
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Fig. 3. Position of meniscus in capillary fit by equation ( 19 1 
with a and h as adjusted parameters. 



from the velocity of the meniscus, w', with the ratio of the 
radius of the circular opening of the gel mold, _R, and that 
of the capillary, r, as u = v'{r/R)'^. The values of R and 
r are measured as 1.745 mm and 0.268 mm, respectively. 
The equilibirium velocity of water flow in the gel is, thus, 
calculated as v^o = 4.4 x 10^^ m/s. The pressure applied 
to the water is Pq = 2.9 x 10'^ N/m^, which corresponds 
to the height of water column of 30 cm, and the thickness 
of the gel is i = 1.0 x lO^'^ m. From these values, the 
friction coefficient of the gel is calculated by equation (pi 
that yields to f — 6.6 x 10^* N s/m*. This agrees with a 
typical value of the friction coefficient of the transparent 
poly (acrylamide) gel evaluated by other experiments |4]. 



3.2 Elastic modulus of gel 

In Figure 4, we show the concentration dependence of the 
Young's modulus and the Poisson's ratio of the gels used 
in our experiment. The least-squares analysis of the results 
yields that the concentration dependence of the Young's 
modulus of the gel is well described by a power law rela- 
tionship with a scaling exponent of 1.7, which is slightly 
smaller than the value expected from the scaling argument 
PI]. The Poisson's ratio is about 0.45 in the dilute gels. It 
decreases with the concentration of the gel and reaches to 
about 0.35 at higher concentration region. The values of 
the Poisson's ratio of the dilute gels are close to those of 
the incompressible materials. On the other hand, the val- 
ues of the Poisson's ratio of the dense gels are rather close 
to those of the metals and the glasses. The Young's mod- 
ulus and the Poisson's ratio of the gel at a concentration 
of 1 M is observed as £; = 1.3 x 10^ N/m^ and a = 0.44, 
respectively. The longitudinal modulus is then calculated 
as K = 4.1 X 10'* N/m^ for the gel at a concentration of 
1 M. Our results are consistent with previously published 
Poisson's ratio of poly(acrylamide) gels and a theory [53J 
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Fig. 4. Concentration dependence of elastic moduli: (a) the 
upper graph shows Young's modulus, E, (b) the lower graph 
shows Poisson's ratio, a. 



The diffusion coefficient of the polymer network is de- 
duced from equation (llj) as D^, = 6.2 x 10^^^ m^/s, which 
agrees with a typical value of the diffusion coefficient of 
the transparent poly(acrylamide) gel evaluated by other 
experiments [51I^E^E5]. 



4 Kinetics of water flow through gel 

When the water ffows in the gel (polymer network) by 
the application of the mechanical pressure, the polymer 
network of gel deforms from the initial position by the 
drag force of the water ffow. We introduce a function, 
u{x,t), that represents the displacement of a point in the 
polymer network from its initial position at x. We assume 
u is always small compared to the characteristic size of 
the gel, L. Under the present definition, u{x, i) = at 
i = 0. We also assume that the macroscopic velocity of 
the water in the polymer network, v{t), and microscopic 
friction coefficient, /, are uniform in space. 



dv 
dx 



0. 



(3) 



Then, we write down the equation of motion of the poly- 
mer network and that of the water in the gel as 



Pgt 



-f{u — v) + K 



dx"^' 



PwV = f{u - v) 



dP 

dx ' 



(4) 
(5) 



where Pg and p„ are the density of the polymer network 
and that of the water, respectively; A dot over the variable 
denotes the time derivative; P{x^ t) represents the pressure 
inside the region that the gel occupied. 

By neglecting the inertia terms in equations (HI) and 
(|5|, we obtain 

d'^u dP ^ 



f{v --«) = - 



dP 

dx ' 



(7) 



A set of equations ([3]) , (|6]) , (|7| corresponds to one dimen- 
sional case under a limit of small polymer volume fraction 
(/)—!■ of a general description for 3 dimensional deforma- 
tion of gels [1S1I17]. 

In our experimental setup of water ffow measurement 
in FigjlJ we see P{x = 0,i) = Pq and P{x = L,t) = 0. 
Taking a spacial average of equation (It]) over the space 
the gel occupied, we obtain 



L 



udx = 



Lf 



(8) 



The time evolution of water ffow in the gel v{t) is deter- 
mined locally from u(x,t) by 



V — ii — Dr 



dx'^ 



(9) 



Here, we use the collective diffusion coefficient of the gel, 
Dc, defined in equation (IT]). The displacement of polymer 
network u{x, t) is, hence, aetermined by solving the equa- 
tion: 



1 



L 



iidx — D, 



d' 



dx'^ 



(10) 



P^ 
Lf 

In the equilibrium state, u{x,t) ^ as i ^ oo. The 
restoring force due to the deformation of the polymer net- 
work is in balance with constant gradient of the pressure. 
In the equilibrium state, therefore, equations (Is]), ([6]), (I?]), 
(|8| give 

d^u{x, t = oo) Po 



dx"^ 



(11) 



DcLf 

The displacement, u(x, t), satisfies the following boundary 
conditions: 



u{L,t) =0, 
9m(0,oo) 



dx 



0. 



(12) 
(13) 



Integration of equation (ITT]) with respect to x with the 
boundary conditons yields 



u{x, oo) = 



Pn 



2D,Lf 



{L^~x^). 



(14) 



Solving equation ( )10| under the above boundary con- 
ditions, we obtain the displacement, u{x,t), as (See Ap- 
pendix) 



(^'*) = Z^^j^[(-l) (^cos— x-1 



exp 



where 



Lf 



(15) 

(16) 
(17) 
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The velocity of the water flow, v(t), is derived from equa- 
tion ^ as 



OO / 9 



l + 2^exp I -— i 

n— 1 



= Woo03(O,e*/^), (18) 



where 03(u, g) is the eUiptic theta function. The position 
of meniscus is expressed as 



OO / 2 

n=l ^ 



where b is the initial position, and 

2 



+ b, 



a 



R 



r 



Vo 



(19) 



(20) 
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5 Discussion 

Now we compare the experimental results with the the- 
oretical predictions by equations (18) and (19 1. The val- 
ues obtained by the experimental measurements are, L = 
1.0 X 10-3 m, K = 4.1 X lO'^ N/m2, and / = 6.6 x lO^^ 
Ns/m^, respectively. 

By applying these values into equation ( 16 1, we obtain 
the characteristic relaxation time t = 1.6 x 10'^ s. The 
equilibrium velocity of water flow in the gel observed by 
the experiments is, as is already given in the section [3J 
Voo = 4.4 X 10~^ m/s. From equation ( 14 1 the deformation 
of gel is estimated as 3.5% which is consistent with our 
assumption of small u. 

The time evolution of the water flow is then calculated 



by equation ( 18 ) with above experimental values. The re- 



sults are given in Figure 5. The normalized velocity of the 
water flow in the gel is also calculated from the slope of 
the flow curve in Figure 3 and shown in Figure 5. We find 
that the time evolution of the velocity of water flow in the 
gel is well explained by the equation (18 1 with a single 
characteristic relaxation time r. 

However, there is a slight discrepancy in the time evo- 
lution curve of the velocity at a short time range (<200s) 
between the theory and the experimental results as one 
can see in figure [3] and [5] The absolute value of the equilib- 
rium velocity determined from the short time data (<200s) 
of the friction experiment is smaller than that determined 
by including the long time data. Figure 5 shows that the 
deference between those two values is roughly 20 %. 

It would be worth noting that the theory predicts that 
the water flow creates a uni-axial deformation of polymer 
network, i.e., the one dimensional concentration gradient 
in the gel, proportional to 



1 



du(x, oo) 
dx 



1 



Pqx 
DcLf 



(21) 



The compression of gel is the cause of the initial rapid 
water flow. It is reported that the friction of the home- 
geneous gel depends on the polymer volume fraction, 0, 



Fig. 5. Water flow velocity in gel: (a) the solid line corresponds 
to the theoretical curve fitting all the data, (b) the dotted line 
corresponds to the theoretical curve expected from the short 
time data (<200s). 



as / oc 0^-^ [Hl21j. The friction, therefore, is expected to 
increase according to the compression of the gel by the 
water flow. The increment of the friction may explain the 
difference between the short time data and long time data. 
The difference between the estimated equilibrium water 
velocities corresponds to 14 % homogeneous compression 
of the sample slab gel, which is larger than the deforma- 
tion which the present theory estimates. Unfortunately, 
we cannot measure the deformation of the gel under the 
water flow with the current experimental setup. 



6 Conclusion 

The kinetics of the water solvent flow through the poly- 
mer network of gel is measured experimentally and de- 
scribed by a simple phenomenological theory. The exper- 
iment shows that the equilibrium state is reached after 
long time and that the friction of the gel shows a large 
value. The time evolution of the velocity of the water flow 
in the gel is calculated on the basis of the collective diffu- 
sion model of the polymer network coupled with the water 
friction, assuming both the water velocity and the friction 
coefficient between the polymer network and the water are 
uniform in the gel. The theory reproduces the time evo- 
lution of the water flow of the experiments for thin dilute 
gel fairly well with a single characteristic relaxation time. 
We are able to deduce the collective diffusion coefficient 
from the water flow experiment on gels. 

However there is a difference between experimental 
value of the equilibrium velocity and the value expected 
from the theoretical calculation by using parameters ob- 
tained from the short time flow measurements and the 
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mechanical measurements. It suggests that the inhomo- 
geneity of the friction and the concentration of the gel un- 
der the water flow might be important for the kinetics of 
the water flow. A better microscopic model is required for 
complete understanding of water flow through gel beyond 
the phenomenological understanding presented in this pa- 
per. 

This work was first suggested by the late Professor Toyoichi 
Tanaka in 1989 during the stay of Y.Y.S. and M.T. in Mas- 
sachusetts Institute of Technology. Y.Y.S. thanks CEA, IPhT 
and Takushoku University, RISE for the financial support. 
M.T. thanks NEDO for the financial support. 



APPENDIX A: Solution of the equation of 
motion 

We define a function W{x, t) as a deviation from the equi- 
librium displacement: 



This indicates i?„ = 0. 



The boundary condition ( 28 ) yields 



dW 
dx 



= Y^[AnCosXnx\{-XlDc)c^^{-XlDJ;) 

n=l 

— w(oo). 



Ao 



This indicates that the constant, A^, should be positive, 
A^ > for n > 1 and Aq = u((X)). Finally, the boundary 
condition ( 29 1 yields 



W{L,t) = ^K(cosA„L + C„A)]exp(-A2D,t)+Aoi + Bo 



n=l 



vdt. 



This indicates CnA = — cosA„i. The following relation- 
ships for Bn, A^, and CnA are obtained: 



u{x, t)- v{t)dt = W{x, t) + u{x, oo). (22) 

We obtain an equation for VF(x, t) from equation (Il0| as 

(23) 



Bn = 0, 

\l > 0, 

CnA = -C0SA„L. 



(30) 
(31) 

(32) 



dt ^ dx"^ 



By assuming W{x^t) — X{x)T{t), the equation (23 1 be- 
comes 

t{t) = -\^D,T{t) 

d^X{x) 



And W{x,t) becomes 

oo 

W{x,t) = J2[McosXnX + CnA)]exp{^XlDJ). (33) 



dx"^ 



= -A"X(x), 



(24) 
(25) 



where A is a constant. The general solution for W{x,t) 
would be 



W{x, t) — \J [A„ cos XnX + Bn siu XnX] 
n=\ 

X exp(-A^i:'ci) + Aqx + Bq. 
The boundary conditions for Wix^i) are 

du _ ^ dW 

dx ^^0 dx 

du ^dW 

dt ^ dt 



Substitution of this solution ( 33 1 into equation ( 23 1 
yields 

oo 

^K(cosA„x-HC„A)](-A2D,)exp(-A2D,i) 
J I ^[A„(cosA„x-f C„a)] 

•'^ n=l 



L 



= 0, 



w(oo). 



u{L,t) = — > W{L,t) 



vdt. 



(26) 

(27) 
(28) 
(29) 



x{~XlDc)eyi^{-XlDct)dx 

OO 

= Dc^ AnXl{- cos A„x) ex^i-XlDct). 



(34) 



The boundary condition ( 27 1 yields 



dW 
dx 



= ^[A„(-A„sinA„a;) -I- B„(A„ cos A„a;)] 
xexp(-A2D,i)|^^^ + Bo 

OQ 

= ^[A„(-A„ X 0) + Bn{Xn X 1)] eM-XlD.t) 

n=l 

+Bo 
= 0. 



The first term of the left handed member (LHMl) is writ- 
ten as follows: 

oo 

LHMl = J2(~^nDc)An{cosXnx)exp{-XlD,t) 

n=l 
oo 

+ Y,i-^lD,){A„CnA)exp{-XlDj). 

The second term of the left handed member (LHM2) is 
expressed as follows: 

LHM2 - - 7 / y2{-XlD,){AnCosXnx) 
X exp{—X^Dct)dx 
-J- Y.^~XlD,){AnCnA)exp{-XlD,t)dx. 
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The first term of LHM2 is expressed as 



L oo 



n=l 



{~XlDc)iAn COS Xnx) exp{-XlDct)dx 



1 , , 1 

= - Y^{-XlDc)^{AnSinXnx)exp{-XlDct) 

1 °° 1 

= - -Y,i-^lDc)T-{AnSmXnL)exp{-XlD,t). 



The second term of LHM2 becomes 



1 
L 

OO 

E 

n=l 



L oo 



n=\ 



{-XID,){A^C^a) eM->^lD,t)dx 



i-XlD,){AnCnA) eM-^lDct). 
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Hence, the equation (34) is reduced to 
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C„A = -cos(7i^) = -(-l)". 
The solution, W{x,t), can be now written as 
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The displacement, u{x,t), is now determined as 

oo 

u(a;,i) = ^A. [(-l)"(cosg„a;-l)]exp(-r„t), (41) 
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where. An, (?«, and Fn are given by 
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